Theory of metastability in simple metal nanowires 
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Thermally induced conductance jumps of metal nanowires are modeled using stochastic Ginzburg- 
Landau field theories. Changes in radius are predicted to occur via the nucleation of surface kinks 
at the wire ends, consistent with recent electron microscopy studies. The activation rate displays 
nontrivial dependence on nanowire length, and undergoes first- or second-order-like transitions as 
a function of length. The activation barriers of the most stable structures are predicted to be 
universal, i.e., independent of the radius of the wire, and proportional to the square root of the 
surface tension. The reduction of the activation barrier under strain is also determined. 

PACS numbers: 05.40.-a, 68.65.La 
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Metal nanowires have attracted considerable interest 
in the past decade due to their remarkable transport and 
structural properties . Long gold and silver nanowires 
were observed to form spontaneously under electron ir- 
radiation 0, @, 0, and appear to be surprisingly sta- 
ble; even the thinnest gold wires, essentially a chain of 
atoms, have lifetimes of the order of seconds at room tem- 
perature 0. Nanowires formed from alkali metals are 
significantly less long-lived, but exhibit striking correla- 
tions between their stability and electrical conductance 
0, 0- That these filamentary structures are stable at 
all is rather counterintuitive lalS li but can be explained 
by electron-shell effects [bl lij. llfij. Nonetheless, these 
nanostructures are only metastable, and understanding 
their lifetimes is of fundamental interest both for their 
potential applications in nanoelectronics and as an inter- 
esting problem in nanoscale nonlinear dynamics. 

In this Letter, we introduce a continuum approach to 
study the lifetimes of monovalent metal nanowires. Our 
starting point is the nanoscale free-electron model 
in which the ionic medium is treated as an incom- 
pressible continuum, and electron-confinement effects are 
treated exactly, or through a semiclassical approxima- 
tion 0, 0, EU- This approach is most appropriate Q 
for studying simple metals, whose properties are deter- 
mined largely by the conduction-band s-electrons, and 
for nanowires of 'intermediate' thickness: thin enough so 
that electron-shell effects dominate the energetics, but 
not so thin that a continuum approach is unjustified. 
The inclusion of thermal fluctuations is modeled using a 
stochastic Ginzburg-Landau classical field theory, which 
provides a self-consistent description of the fluctuation- 
induced thinning/growth of nanowires. 

Our theory provides quantitative estimates of the life- 
times for alkali nanowires with electrical conductance G 
in the range 3 ^ G/Go ^ 100, where Gq = 2e 2 /h is the 
conductance quantum. In addition, we predict a univer- 
sality of the typical escape barrier for a given metal, in- 
dependent of the wire radius, with a value proportional 
to %fa, where a is the surface tension of the material. 
Our model can therefore account qualitatively for the 



large difference in the observed stability of alkali vs. no- 
ble metal nanowires. It also predicts a sharp decrease of 
the escape barrier under strain. 

We consider a cylindrical wire suspended between two 
metallic electrodes, with which it can exchange both 
atoms and electrons. The resulting energetics is de- 
scribed through an ionic grand canonical potential 



a, 



(1) 



where fl e is the free energy for a fixed number Af a of 
atoms and fj, a is the chemical potential for a surface atom 
in the electrodes. In the Born-Oppenheimer approxima- 
tion, fl e is just the electronic grand-canonical potential, 
and can be written as a Weyl expansion plus an electron- 
shell correction 



-UjV + crS - 



dzV(R(z)), 



(2) 



where V, S, and L are, respectively, the volume, sur- 
face area, and length of the wire, ui and a are material- 
dependent coefficients, and V(R), shown in Fig. ^ is 
a mesoscopic electron-shell potential .1QJ that describes 
electronic quantum-size effects. 
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FIG. 1: Electron-shell potential V(R) at zero temperature. 
The top axis shows the conductance values of the most stable 
wires in units of the conductance quantum, Go = 2e 2 jh. 
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In the presence of thermal noise, a wire's radius will 
fluctuate as a function of time t and position z measured 
along the wire's axis: R(z, t) = R + <f>(z, t) for a wire of 
radius R at zero temperature. The wire energy Q may 
be expanded as 



n a [R,cf)} = n a (R) + H[<f>] ) 



(3) 



where H.[(f)] is the energy of the fluctuations. Keeping 
only the lowest-order terms in d z <p, one finds 



n[tf>] 



dz 



(4) 



where k = 2iraR and U(4>) is an effective potential. 

A (meta)stable nanowire is in a state of diffusive equi- 
librium: 

-v. (5-^^), 0) 

where /x cy i is the chemical potential of a surface atom in 
a cylindrical wire and V a is the volume of an atom. 
Using Eq. JSJ) in Eq. Q), one finds the effective potential 



u(q>) = v(R+<j>)-v(R)-^- 



£) (6) 
2R dR V ' 



A stable wire must satisfy U'{0) = and U"(0) > 0. 
The first condition is satisfied automatically. The second 
condition is equivalent to the requirement dfi cy \/dR > 0, 
and was previously used to determine the linear stability 
of metal nanowires [|J . The most stable wires correspond 
to the minima of V(R) (cf. Fig. QJ; however, the stable 
zones span finite intervals of radius about the minima . 

The radius fluctuations <fr(z, t) due to thermal noise can 
be treated as a classical field on [0, L) , with dynamics gov- 
erned by the stochastic Ginzburg-Landau (GL) equation 



d<Kz,t) 
dt 



d 2 

'dz 



dU 



+ (2Ty/ 2 az,t), 



(7) 



where £(z, i) is unit-strength spatiotemporal white noise. 
The zero-noise dynamics is gradient, that is, = —8TL/5<j) 
at zero temperature. In J7J, time is measured in units of 
a microscopic timescale describing the short- wavelength 
cutoff of the surface dynamics 0, |?J which is of order the 
inverse Debye frequency . 

At nonzero temperature, thermal fluctuations can 
drive a nanowire to escape from the metastable configu- 
ration (f> = 0, leading to a finite lifetime of such a nanos- 
tructure. The escape process occurs via nucleation of a 
"droplet" of one stable configuration in the background 
of the other, subsequently quickly spreading to fill the 
entire spatial domain. A transition from one metastable 
state to another must proceed via a pathway of states, 
accessed through random thermal fluctuations, that first 



goes "uphill" in energy from the starting configuration. 
Because these fluctuations are exponentially suppressed 
as their energy increases, there is at low temperature a 
preferred transition configuration (saddle) that lies be- 
tween adjacent minima. The activation rate is given in 
the T — > limit by the Kramers formula 



T - r exp(-A£/T) 



(8) 



Here AE is the activation barrier, the difference in energy 
between the saddle and the starting metastable configu- 
ration, and Tq is the rate prefactor. 

The quantities AE and Tq depend on the microscopic 
parameters of the nanowire through k and the details of 
the potential ©, on the length L of the wire, and on 
the choice of boundary conditions at the endpoints z — 
and z — L. Simulations of the structural dynamics un- 
der surface self-diffusion 01 suggest that the connection 
of the wire to the electrodes is best described by Neu- 
mann boundary conditions, d z 4>\o,L = 0. These bound- 
ary conditions force nucleation to begin at the endpoints, 
consistent with experimental observations |J|. 

The saddle configurations are time-independent solu- 
tions of the zero- noise GL equation 0] , and can be ob- 
tained by numerical integration of Eq. Q at T = 0. 
However, we find that for many of the metastable wires, 
the effective potential U (4>) can be approximated locally 
by a cubic potential 



C/W(0) = - a ± + |(0 ± ) 3 , 



(9) 



where 



=F 4> and a, f3 > 0. The potential 



(U^) biases fluctuations toward smaller (larger) 
radii. It is useful to scale out the various constants 
in the model by introducing the dimensionless variables 
x = z/Lq and u = (/3/a) 1 / 2 ^, where L = k 1 / 2 /(a/3) 1 / 4 
and Eq = are characteristic length and 

energy scales. The energy functional then becomes 



H[u] 
E 



~ 2 {u) 



dx. 



(10) 



where £ = L/Lq and u' — du/dx. 

Metastable and saddle configurations are stationary 
functions of Eq. QlOjl. and therefore obey the Euler- 
Lagrange equation u" = — 1 + u 2 . With Neumann bound- 
ary conditions, the uniform stable state is the constant 
state u s = +1, and there exists a uniform unstable state 
u u = —I. We will see that the latter is the saddle for 
£ < t c = Tr/y/2. At £ c a transition occurs 13], and above 
it the saddle is nonuniform. It consists of an 'instanton' 
localized at one end of the wire, and is given by [l4| 



u{x) 



dn 2 



(ii) 
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where dn(- | m) is the Jacobi elliptic dn function with pa- 
rameter m, with to 1, and £(m) = to 2 — m + 1. Its 
half-period is given by K(m), the complete elliptic inte- 
gral of the first kind, a monotonically increasing function 
of to. Eq. satisfies the Neumann boundary condition 
when 



quency I'd)- For £ < 



£ = V2£(to) 1/4 K(to), 



(12) 



which (taking to — > 0+) leads to 4 = tt/V2. As £ -> £+, 
dn(x|0) — > 1, and the nonuniform saddle reduces to the 
uniform state u u = —1. This is the saddle for all £ < £ c . 

When the saddle is constant (£ ^ £ c ), the activation 
barrier scales linearly with (reduced) length £: AE/Eq = 
(4/3)£. Above l c , it is expressible in terms of the com- 
plete elliptic integrals of the first kind K(to) and the 
second kind E(to): 



AE 

E{) 
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2 
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(2 — m)(l — to) 



K(to) 



(13) 



As £ — > oo (corresponding to to — > 1~), AE/Eq — > 
12-\/2/5. More generally we denote the asymptotic value 
linif^oo AE(L) = AEoo. The activation barrier for the 
entire range of £ is shown in Fig. 

Calculation of the prefactor To in the Kramers tran- 
sition rate formula is a much more involved matter. It 
generally requires an analysis of the transverse fluctua- 
tions about the extremal solutions. The general method 
for determining T has been discussed elsewhere 0, ^| ; 
here, we just present results (in units of the Debye fre- 
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FIG. 2: The activation energy AE as a function of the wire 
length L, for the cubic potential with Neumann boundary 
conditions (top). The dashed line indicates the critical wire 
length L c at which the transition takes place. The bottom 
panel shows the prefactor To, and the inset displays the ac- 
tivation barrier for the full potential U((j>) for kpRo = 12.79, 
exhibiting a succession of first order phase transitions. 
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we find 

1 sinh(£ v / 2) 



tt sin 



(14) 



which diverges as £ — > £ c ~ , with a critical exponent of 
1/2. The divergence arises from a soft mode; one of the 
eigenmodes corresponding to small fluctuations about 
the saddle has vanishing eigenvalue at £ c - This diver- 
gence, and its meaning, are discussed in detail in [T^ . 
For £ > £ c , the prefactor is 



2 - to + 2\/4to 2 - 



1 



4tt^(to) 3 / 8 



(1 — to) sinh 


2£(to) 1 /4K(to) 




£(m)E(m) - |(1 - m)(2 - m)K(m) 



(15) 



This also exhibits a divergence with a critical exponent 
of 1/2 as £ -> £ c + . The prefactor over the entire range of 
£ is shown in Fig. [21 

The second-order-like transition in activation behavior 
exhibited in Figgis interesting, but generally holds only 
for transitions where the potential U(<f>) can be locally 
approximated by a smooth potential of quartic or lower 
order .14|. For some of the minima of Fig. ^ this is 
not the case, and the wire instead exhibits one or more 
first-order- like transitions 0], as shown in the inset of 
Fig.0 

Figure shows the activation barrier AEoo as a func- 
tion of radius R for a typical metastable wire, correspond- 
ing to the conductance plateau at G = 17 Go in Au. 
Very good agreement is found between the numerical re- 
sult for the full potential © (solid curve) and the result 
from Eq. 1)13(1 using the best-fit cubic polynomial 

F [nN] 
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FIG. 3: The activation energy AEoo as a function of radius 
for a typical stable zone in Au. Solid curve: numerical result 
for the full potential U(<f>), Eq. ©; dashed curve: result from 
Eq. 1131 using the best cubic-polynomial fit to U (4>). The 
wire radius is related to the tensile stress (upper axis). 
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TABLE I: The lifetime r (in seconds) for various cylindrical 
sodium nanowires at temperatures from 75K to 125K. Here 
G is the electrical conductance of the wire, L c is the criti- 
cal length above which the lifetime may be approximated by 
r ~ v^ 1 exp(A_Boo/T), and A_Eoo is the activation energy for 
an infinitely long wire. Note that wires shorter than L c are 
predicted to have shorter lifetimes. 



G 

[Go] 


L e 

[A] 


[meV] 


r[s] 


75 K 


100 K 


125 K 


3 


2.8 


250 


4 x 10 5 


2 


5 x 10" 3 


6 


4.3 


200 


7 


3 x 10~ 3 


3 x 10~ 5 


17 


5.0 


260 


7 x 10 5 


3 


8 x 10~ 3 


23 


6.1 


230 


2 x 10 3 


0.2 


9 x 10~ 4 


42 


7.2 


250 


2 x 10 5 


1 


io- 3 


51 


6.8 


190 


1 


8 x 10" 4 


10" 4 


67 


18.8 


180 


0.6 


5 x 10 -4 


7 x 10" 6 


96 


11.4 


250 


10 5 


0.8 


3 x 10" 3 



(dashed curve). Under strain, R varies elastically; the 
corresponding stress in the wire is shown on the upper 
axis. A stress of a fraction of a nanonewton can signifi- 
cantly change the activation barrier, and even change the 
direction of escape. The maximum value of AE^ occurs 
at the cusp, where the activation barriers for thinning 
and growth are equal. 

The most stable structures, corresponding to the max- 
imum values of AEqc, occur at (or near) the minima 
of the electron-shell potential, V(R) (Fig. 0. The life- 
times of these equilibrated structures are limited by thin- 
ning, since the total energy of the wire is lowered by 
reducing its volume. We thus fit the effective poten- 
tial at these minima to the form ZJ( - ). Table [Q lists 
critical lengths L c , activation barriers AE^, and life- 
times t = 1/r, Eq. (JSJl, at various temperatures for Na 
nanowires. (Only the minima that are well-fit by [A - ) 
are shown.) The temperature dependence of r shows that 
the lifetime of Na nanowires drops below the threshold for 
observation in break-junction experiments as the temper- 
ature is increased from 75K to 125K. This behavior can 
explain the observed temperature dependence of conduc- 
tance histograms for Na nanowires which show clear 
peaks at conductances near the predicted values at tem- 
peratures below 100K, but were not reported at higher 
temperatures. The increase of L c with G, shown in Ta- 
ble[H may also explain the observed exponential decrease 
in the heights of the conductance peaks with increasing 
conductance , since the thicker contacts are more likely 
to be shorter than L c , and hence to have exponentially 
reduced lifetimes. 

An important prediction given in Table [Q is that the 
lifetimes of the most stable nanowires, while they do ex- 



TABLE II: Fermi energy, surface tension, and typical criti- 
cal length and activation barrier for various alkali and noble 
metals. The surface tension values are extrapolations to zero 
temperature from Ref. [l6||. 



Metal 


Li 


Na 


K 


Rb 


Cs 


Cu 


Ag 


Au 


e F [eV] 
a [N/m] 


4.74 
0.52 


3.24 
0.26 


2.12 
0.14 


1.85 
0.12 


1.59 
0.09 


7.00 
1.78 


5.49 
1.24 


5.53 
1.50 


L c /R 
AE X [meV] 


0.67 
290 


0.71 
200 


0.81 
150 


0.84 
140 


0.88 
120 


0.83 
530 


0.88 
440 


0.97 
490 



hibit significant variations from one conductance plateau 
to another, do not vary systematically as a function of 
radius; the activation barriers in Table [I] vary by only 
about 30% from one plateau to another, and the wire 
with a conductance of 96Go has essentially the same life- 
time as that with a conductance of 3Gq. In this sense, the 
activation barrier is found to be universal: in any con- 
ductance interval, there are very short-lived wires (not 
shown in Table P) with very small activation barriers, 
while the longest-lived wires have activation barriers of a 
universal size 

AEod ^ o.6 — , (16) 

V m e J 

depending only on the surface tension of the material. 
Here m e is the conduction-band effective mass, which 
is comparable to the free-electron rest mass. The fact 
that the typical activation energy l|16fl is independent 
of R is a consequence of the virial theorem: Since the 
instanton is a stationary state of Eq. Q, the bending 
energy (^(9 z </>) 2 ) is proportional to (U(cf))). Since k ~ 
o~ R and V ~ 1 /R 01 > this implies that the characteristic 
size of the instanton L c ~ \faR and AEoo ~ ^/a. 

Table [H] lists typical activation barriers and critical 
lengths for various alkali and noble metals. It shows that 
noble metal nanowires should have much longer lifetimes 
than alkali metal nanowires, due to their larger surface 
tension coefficients. This prediction is consistent with 
experimental observations 0, 0, [3 IE IE > although our 
estimated activation barriers for noble metal nanowires 
are still too small to account for their observed stability 
at room temperature. This discrepancy may stem from 
the neglect of d-electrons in our model (except inasmuch 
as they enhance a compared to the free-electron value), 
or due to the presence of impurities which passivate the 
surface, thereby raising the activation barrier above its 
intrinsic value. 

J.B. and C.A.S. acknowledge support from NSF Grant 
No. DMR0312028. D.L.S. acknowledges support from 
NSF Grant Nos. PHY0099484 and PHY0351964. 



■5 



[1] N. Agra'ft, A. Levy Yeyati, and J. M. van Ruitenbeek, 

Phys. Rep. 377, 81 (2003). 
[2] Y. Kondo and K. Takayanagi, Phys. Rev. Lett. 79, 3455 

(1997). 

[3] V. Rodrigues, J. Bettini, A. R. Rocha, L. G. C. Rego, 

and D. Ugarte, Phys. Rev. B 65, 153402 (2002). 
[4] Y. Oshima, Y. Kondo, and K. Takayanagi, J. Electron 

Microsc. 52, 49 (2003). 
[5] R. H. M. Smit, C. Untiedt, and J. M. van Ruitenbeek, 

Nanotech. 15, S472 (2004). 
[6] A. I. Yanson, I. K. Yanson, and J. M. van Ruitenbeek, 

Nature 400, 144 (1999). 
[7] A. I. Yanson, J. M. van Ruitenbeek, and I. K. Yanson, 

Low Temp. Phys. 27, 807 (2001). 



[8] F. Kassubek, C. A. Stafford, H. Grabert, and R. E. Gold- 
stein, Nonlinearity 14, 167 (2001). 
[9] C.-H. Zhang, F. Kassubek, and C. A. Stafford, Phys. 
Rev. B 68, 165414 (2003). 

[10] J. Biirki, R. E. Goldstein, and C. A. Stafford, Phys. Rev. 
Lett. 91, 254501 (2003). 

[11] C. A. Stafford, D. Baeriswyl, and J. Biirki, Phys. Rev. 
Lett. 79, 2863 (1997). 

[12] P. Hanggi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 
62, 251 (1990). 

[13] R. S. Maier and D. L. Stein, Phys. Rev. Lett. 87, 270601 
(2001). 

[14] D. L. Stein, J. Stat. Phys. 114, 1537 (2004). 

[15] E. M. Chudnovsky, Phys. Rev. A 46, 8011 (1992). 

[16] W. R. Tyson and W. A. Miller, Surf. Sci. 62, 267 (1977). 



